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Abstract 

The integrability problem consists in finding the class of functions a first in- 
tegral of a given planar polynomial differential system must belong to. We recall 
the characterization of systems which admit an elementary or Liouvillian first 
integral. We define Weierstrass integrability and we determine which Weierstrass 
integrable systems are Liouvillian integrable. Inside this new class of integrable 
systems there are non-Liouvillian integrable systems. 
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1 Introduction 

It is not always possible and sometimes not even advantageous to explicitly write the 
solutions of a system of differential equations in terms of elementary functions. In fact, 
Poincare begun the qualitative theory of differential equations to well-understand the 
behavior of the solutions of a differential system without their explicit knowledge. For 
Poincare, it is thus necessary to study the functions defined by the differential equations 
by themselves and without bringing them back to simpler functions. These thoughts 
induced Poincare to tackle the study of differential equations beyond an essentially dif- 
ferent point of view from his predecessors. His study provokes a conceptual change on 
the understanding of differential equations. Sometimes, though, it is possible to find 
elementary functions that are constants on solution curves, that is, elementary first 
integrals. These first integrals allow to occasionally deduce properties that an explicit 
solution would not necessarily reveal, see for instance [T2]. This thought originated 
the modern integrability theory of differential equations that tries to respond to the 
natural question: When does a system of differential equations have a first integral that 
can be expressed in terms of "known functions" an how does one find such an integral? 
The answer when the "known functions" are the elementary functions (i.e. functions 
expressible in terms of exponentials, logarithms and algebraic functions) was given in 
|12j . and when the "known functions" are the Liouvillian functions (i.e. functions that 
are built up from rational functions using exponentiation, integration, and algebraic 
functions) was given in 13J. In these two cases it is given the form of an integrating 
factor if the system has these type of first integrals. In this paper we extend the results 
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presented in [T2] and [13]. 



In elementary courses on differential equations we consider systems of the form 

dx „ , dv „ , , 

i = — = P(as, y), y = = Q(z, y), (1.1) 

where P and Q are polynomials in C[a;,?/], C being the complex numbers. Through- 
out this paper we will denote by m - maxjdegP, deg Q} the degree of system (|1.1[) . 
Obviously, we can also express system (jl.lj) as the differential equation 

dy = Q(x,y) 

dx P{x,y) ' [ ' 1 

We learn that although we cannot always explicitly solve this system, we are oc- 
casionally able to find first integrals, that is, nonconstant functions H(x,y), ana- 
lytic on some nonempty open set in C 2 , that are constant on the solution curves 
in this set. To do this we consider the differential form Q(x,y)dx — P{x,y)dy = 0. 
If dP/dx = —dQ/dy, then H(x,y) = j Qdx — Pdy will be a first integral. If 
dP/dx ^ —dQ/dy, we are taught ad hoc methods to find an integrating factor, 
that is a function R(x,y) such that d(RP)/dx = —d(RQ)/dy. In case we can find 
such function R, H(x, y) = J RQdx — RPdy will be a first integral. For example, if 
(dQ/dx + dP/dy)/P is independent of y, then R = exp(J (dQ /dx + dP/dy)/ Pdx) will 
be an integrating factor. 



2 Integrability problem 

We recall that the integrability problem consists in finding the class of functions a first 
integral of a given system (|1.1[) must belong to, see [2]. For instance in [11], Poincare 
stated the problem of determining when a system (| 1 . ljl has a rational first integral. 
The works of |12j and [13] go in this direction since they give a characterization of when 
a polynomial system (|1.1|) has an elementary or a Liouvillian first integral. A precise 
definition of these classes of functions is given in [T2l [T3] . An important fact of their 
results is that invariant algebraic curves and exponential factors play a distinguished 
role in this characterization. Moreover, this characterization is expressed in terms of 
the inverse integrating factor. Now, we state some results related to integration of a 
system (|l.lj) by means of elementary and Liouvillian functions. 

Theorem 2.1 [12j // system has an elementary first integral, then there exists 

wo, oj±, . . . ,u> n algebraic over the field C(x, y) and ci, C2, . . . , c„ in C such that the ele- 
mentary function 

n 

H = oj Q + ^2c l \n(w l ), (2.3) 

i=l 

is a first integral of system il.l]) . 

The existence of an elementary first integral is intimately related to the existence 
of an algebraic inverse integrating factor, as the following result shows. 

Theorem 2.2 |12j // system (1.1)) has an elementary first integral, then there is an 
inverse integrating factor of the form 

( A(x,y) \ 1/N 
\B(x,y)J ' 
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where A, B 6 C[x, y] and N is a nonnegative integer number. 



In the work [3], the systems with a (generalized) Darboux first integral, that 
is, with a first integral of the form 

where /i, g^, /ii £ C[x, y], Xi, \i% £ C for V« and £ N for i = 1, . . . , £, are studied and 
the following result is accomplished. 

Theorem 2.3 3J If system hl.l]) has a (generalized) Darboux first integral of the form 
[2-4% then there is a rational inverse integrating factor, that is, an inverse integrating 
factor of the form 

A{x,y) 
B(x,yY 

where A, B e C[x, y]. 

Unfortunately, not all the elementary functions of the form (|2 .3|) are of (generalized) 
Darboux type. That's why we can find systems with an elementary first integral and 
without a rational inverse integrating factor. We remark that both Theorems 12.21 and 
12.31 give necessary conditions to have an elementary or (generalized) Darboux, respec- 
tively, first integral. The reciprocal to the statement of Theorem l2.2l is not necessarily 
true. But the reciprocal to the statement of Theorem 12.31 is true as we will see later. 



The following Theorem 12.41 ensures that given a (generalized) Darboux inverse in- 
tegrating factor, there is a Liouvillian first integral. The Liouvillian class of functions 
contains the rational, Darboux and elementary classes of functions. Singer gives in [13] 
the characterization of the existence of a Liouvillian first integral for a system (|l.ip by 
means of an integrating factor. 

Theorem 2.4 [13j System il.l]) has a Liouvillian first integral if, and only if, there 
is an inverse integrating factor of the form V = exp j J*^'^ v\> where n is a rational 
1-form such that dr\ = 0. 

Taking into account Theorem 12.41 Christopher in [3] gives the following result, 
which makes precise the form of the inverse integrating factor. 

Theorem 2.5 [5j If the system il.l]) has an inverse integrating factor of the form 
V = exp | I(xoy ) w h ere V * s a rational 1-form such that drj = 0, then there exists 
an inverse integrating factor of system 1 1.1]) of the form 

V = exp{D/E}l[cl% (2.5) 

where D, E, and the Ci are polynomials in x and y and Zj G C. 

We notice that d — and E = are invariant algebraic curves and exp{D/E} is 
an exponential factor for system see for instance [3]. Theorem 12 .51 states that the 

search of Liouvillian first integrals can be reduced to the search of invariant algebraic 
curves and exponential factors. A result to clarify the easiest functional class of the first 
integral once we know the inverse integrating factor is a straightforward consequence 
of Theorem 12.51 and is the reciprocal of Theorem [ 
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Corollary 2.6 If system has a rational inverse integrating factor, then the sys- 

tem has a (generalized) Darboux first integral. 

The proof is based in that if V is a rational inverse integrating factor of system 
p.ip then, 77 = (Q(x,y)dx — P(x,y)dy)/V(x,y) is a rational 1-form such that dn = 0. 

Since H = exp | J^^yo) ^} ^ s a ^ rs ^ m tegral of system by Theorem 12.51 H is a 

(generalized) Darboux hrst integral. 

In [TO], Painleve proved the following result, see also [3(7] and references therein. 

Theorem 2.7 10] A differential system il.ll) has a first integral of the form 

I(x, y) = (y- 9l (x)r (y ~ 92(x)) a2 • . . (y - gt(x)) a *h(x) , (2.6) 

where gj(x) are unknown particular solutions of 11.2]) , h(x) is an unknown function of 
x and the on are unknown constants such that Y\i=i a i 7^ 0> if and only if it has an 
integrating factor of the form 

„ ,/ x a(x)S(x, y) .„ 

M(x, y) = V ; V ' y) — - , (2.7) 

(y - 9i{x)){y ~ 92{x)) ...{y- g t {x)) 

where S{x, y) is polynomial in the variable y of degree t — m — 1. Moreover, 

(a) if the system has two different integrating factors Mi and Mi of the form \2. 7| ) 
with M2 /M\ nonconstant, then there exists a change of variable that is rational 
in the variable y which transforms the equation U.2\) into a Riccati equation. 



(b) if the differential system has only one integrating factor of the form {2.7}) , then 
the particular solutions gi{x) from the ansatz \2. 6|j are calculated algebraically 
and h(x) is given by a logarithmic quadrature. 

As usual we define C[[x]] the set of the formal power series in the variable x with 
coefficients in C and C[y] the set of the polynomials in the variable y with coefficients 
in C. A polynomial of the form 



n 

£« 

i=0 



is called a formal Weierstrass polynomial in y of degree n if and only if a n {x) = 1 and 
cti(O) = for i < n. A formal Weierstrass polynomial whose coefficients are convergent 
is called Weierstrass polynomial, see [I] . In a natural generalization we call Weierstrass 
rational function a function which is a quotient of sums of Weierstrass polynomials. 
We say that a system (ll.l|) is Weierstrass integrable if system (jl.l[) admits an inverse 
integrating factor of the form 

V = exp{D/E}]Jcl i , (2.8) 

where D, E, and the (7, are Weierstrass polynomials and U G C. In this sense, the 
systems which are Liouvillian integrable are included in the set of systems which are 
Weierstrass integrable because they have an inverse integrating factor of the form (|2.8[) . 
However, there are systems which are Weierstrass integrable which are not Liouvillian 
integrable, see Example 2 in [5]. The systems with a first integral of the form (|2.6|) are 
Weierstrass integrable because they have an integrating factor of the form (12. 7|) . 

The following theorem determines some Weierstrass integrable systems which are 
Liouvillian integrable. 
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Theorem 2.8 If a differential system (1.1)) has a first integral of the form \2.6\) , and 
at least one algebraic solution, then it has a Liovillian first integral and therefore a 
(generalized) Darboux inverse integrating factor of the form &2.5\) . 



Proof. If the system has a first integral of the form (|2.6[) , then by Theorem 12.71 either 
there exists a rational Weierstrass change which transforms the system of differential 
equations (jl.l)) into a Riccati equation (see statement (a)) or all the particular solu- 
tions gi(x) from the ansatz (|2.6[) are calculated algebraically and h(x) is given by a 
logarithmic quadrature (see statement (b)). In the first case, as we have an algebraic 
solution, the Riccati equation can be solved by quadratures and the system has a Li- 
ouvillian first integral. In the second case all the curves y — gi{x) = are algebraic 
curves and h(x) is given by a logarithmic quadrature, which implies that the inverse 
integrating factor (|2.7[) either is a rational integrating factor if a(x) is a polynomial or 
is a (generalized) Darboux integrating factor if a{x) is not a polynomial. In both cases 
by Corollary |2.6l or Theorem l2.5l the system has a Liouvillian first integral. | 

In the following examples we will see the existence of planar polynomial systems 
which are not Weierstrass integrable. 

In [5j it is proved that all the rational Abel differential equation known as solvable in 
the literature can be reduced to a Riccati differential equation or to a first-order linear 
differential equation through a change with a rational map. Several examples of Abel 
differential equations which are not Weierstrass integrable appear in the Appendix of 
[S]. For instance the Abel differential equation 



^=y 3 -2xy 2 , (2.9) 
ax 



has the following first integral 
H(x,y) = 



xAi (x 2 - + Ai (l,x 2 
xBi (x 2 - ±) +Bi(l,; 



•2 _ I 

V 



where Ai and Bi are a pair of linearly independent solutions of the Airy differential 
equation u" — zlo. However, the change X — x 2 — 1/y, Y — x transforms equation 
(|2.9p into the Riccati equation dY/dX — Y 2 ~X and in these new variables the system 
is Weierstrass integrable. The same happens for all the other cases given in [5]. 

In [9] it is presented a new algorithm to detect analytic integrability or a singular 
expansion of the first integral around a singular point for planar vector fields. It is also 
studied the following example 

x = -y, y = ax + by + y 2 , (2.10) 

and it is proved that it has a first integral of the form 

oo x 

where Pk is a polynomial of degree < k. Hence, system (|2.10|) is not Weierstrass 
integrable and for this case it is unknown if there exists a change of variables which 
transforms the system into a Weierstrass integrable one. 
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